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Abstract. Spontaneous breaking of chiral symmetry in QCD has traditionally been inferred 
indirectly through low-energy theorems and comparison with experiments. Thanks to the 
(— I ' understanding of an unexpected connection between chiral Random Matrix Theory and chiral 

O |. Perturbation Theory, the spontaneous breaking of chiral symmetry in QCD can now be 

shown unequivocally from first principles and lattice simulations. In these lectures I give an 
introduction to the subject, starting with an elementary discussion of spontaneous breaking of 
^ ( — i global symmetries. 



1. Introduction 



> 

in 

£\\ ' Many tend to think of spontaneous breaking of global symmetries in quantum field theory as 

something rather simple, almost trivial. To illustrate that this is certainly not so, let us first 
briefly review the standard textbook treatment of this phenomenon. Typically, one considers 
a field theory of a real scalar field with a potential that, appropriately for this school, can be 
dubbed a "sombrero potential" , 

^ . VQ) = - \m 2 ^x) 2 + ^{xf . (1) 

H ■ 

From a perspective of classical physics this looks unstable. Certainly a constant field 
configuration of <j){x) = cannot possibly be a good starting point for a perturbative expansion. 
Instead, one considers the two minima of the potential at 

4> = ±v = ± \ -m (2) 
V A 

and says that a consistent treatment of such a field theory must be based on one vacuum that 
is either at <j) = v or at <j) = —v, which one being undecidable. Here a comparison is often made 
to a classical picture, like a pencil standing exactly on its tip: such a system is rotationally 
symmetric, but any small perturbation will evidently make the pencil fall in the direction of the 
perturbation. Spontaneous symmetry breaking! Indeed, the pencil will now lie flat on the table, 
and it will have picked just one random direction, the one that was induced by the perturbation, 
a perturbation we can make as small as we like. What about the potential energy of the pencil? 
It was first converted to kinetic energy, and then, on impact with the table, to heat. Such a 



system is dissipative: If it were not, the pencil would simply bounce back up. There may issues 
such as accuracy of initial conditions and so on, but this last point illustrates the difficulty with 
a simple picture of spontaneous symmetry breaking. 

To investigate this a little more closely, let us proceed to quantum mechanics and glue two 
harmonic oscillator potentials together, separated by a distance 2a0. The potential is thus 

V(x) = imw 2 (\x\ - a) 2 (3) 

The ground state of this quantum mechanical system must have no nodes, and we can almost 
trivially identify the correct ground state by drawing two Gaussian bumps centered at x = ±a, 



i/j± (x) ~ exp 



1 / ^2 
--mu{x =p a) 



(4) 



which are, individually, ground states of either of the two harmonic oscillators. If we smoothly 
join the two at x = so that we construct, roughly, 

ipo(x) ~ ip+(x) +ip-(x) (5) 

this will be a very good approximation to the true ground state. What about the 1st excited 
state? It should have one node, and again we can easily guess its form. Suppose we glue the 
two wave functions ip±{x) together, but with a twist: we flip the sign of one of them to get 

i>l(x) ~ ip+(x) - ip-(x) (6) 

This state will have a slightly higher energy than ipo(x), but it is clear from this simple 
construction that ipo(x) and ipi(x) are almost degenerate. The higher excited states can be 
visualized similarly, by gluing two wave functions together around x = 0, more and more excited 
states, and alternatingly adding and subtracting combinations. 

Does such a system display spontaneous symmetry breaking? It is well known that in 
this quantum mechanical case one of the two apparently distinct vacua will never be selected. 
Although the first excited state is almost degenerate with the ground state, there will always 
be a finite energy difference E\ — Eq. If we prepare a state that initially is localized in only 
one of the two wells, it will tunnel to the other side. The tunneling time goes like (E% — Eq) . 
Only when the two minima are infinitely far from each other does that tunneling time go to 
infinity. So our naive picture of spontaneous symmetry breaking can clearly not be valid for 
such a simple quantum mechanical system of just one degree of freedom. One cannot contain 
the probability density on just one side of x = 0. 

These two examples, one classical and one quantum mechanical, illustrate that the 
phenomenon of spontaneous symmetry breaking in quantum field theory cannot be quite as 
simple as the intuition based on the sombrero potential seems to suggest. Nevertheless, as is so 
often the case in physics: the argument may not be quite right, but the answer is correct. There 
is no spontaneous breaking of symmetry in quantum mechanical systems of a finite number of 
degrees of freedom. But there can be spontaneous breaking of symmetries in quantum field 
theory. The distinction here comes precisely from the infinite number of degrees of freedom, 
the infinite number of qauntum field theory modes. Although each of them individually can 
prevent spontaneous breaking of symmetry by tunneling, this mechanism can be blocked when 
an infinite number of modes have to tunnel. Subtleties in the arguments have been discussed on 
and off in the literature, see for instance refs. [21 E]- For a more recent discussion, with plenty 
of references to earlier literature, see ref. [3]. 



1 This is also a standard textbook example, see e.g., ref. [I] for a nice and clear discussion. 



To make the confusion complete, let us also point out another related phenomenon in quantum 
field theory: tunneling from vacuum to another. We have just argued that in quantum field 
theory the tunneling between degenerate vacua can be prohibited due to the fact that an infinite 
number of degrees of freedom have to simultaneously align. Yet, if we tilt the sombrero potential 
slightly so that one vacuum really becomes energetically favored it is also common knowledge 
that in quantum field theory we do have tunneling to the true vacuum. How can this be? It 
turns out to be another example of a superficial argument that is basically correct nevetheless, 
and there is an endless amount of literature on this subject. Suffice it to say here that when 
there is a genuine non-degeneracy of vacua, one can have bubble formation and a first order 
phase transition where the true vacuum gobbles up the false vacuum. For a nice discussion of 
the essential physics of this situation see, e.g., ref. [SJ. 

Of particular importance is the phenomenon of spontaneous breaking of continuous global 
symmetries because it is linked to what is known as Goldstone's Theorem. This theorem 
implies the appearance of one massless mode associated with each broken generator of the 
continuous group of symmetries. Again there are caveats, and there are several. One is space- 
time dimensionality: a massless scalar in two space-time dimensions is not well defined in an 
infinite volume due to infrared singularities, and indeed there is no Goldstone phenomenon in 
two space-time dimensions [6j. In higher dimensions the theorem is incredibly strong, as it says 
that there is no mass gap, an otherwise forbiddingly difficult question to tackle. The Goldstone 
Theorem is most concisely stated in the operator language. Since the symmetry is continuous, 
let us consider the associated conserved current operator. Ordinarily the vacuum |0) is assumed 
to be annihilated by the corresponding conserved charge operator Q (the vacuum is not charged 
under Q). If the vacuum is not annihilated by Q, the symmetry is spontaneously broken. 
Remarkably, this implies the existence of an associated massless mode. This is the content of 
Goldstone's Theorem. 

In QCD, the spontaneous breaking of chiral symmetry is a profound phenomenon that has 
been inferred indirectly over a long period of both theoretical and experimental development, 
starting from the days when the fundamental theory of strong interactions was not even known. 
In hindsight, it is incredible that it could be gotten at without a correct understanding of the 
underlying microscopic mechanism. Today we're much better off: we know the fundamental 
Lagrangian of QCD and it is "just" a matter of checking whether chiral symmetries are broken 
spontaneously or not. In practice, this is not simple at all. In fact, the question of showing 
spontaneous breaking of chiral symmetry in QCD and QCD-like theories is a challenge at the 
level of showing confinement and the existence of a mass gap in pure Yang-Mills theory. It is 
one of the remarkable successes of the numerical non-perturbative approach to QCD - lattice 
QCD - that this now has been shown beyond any doubt. An essential ingredient in this comes 
from the new understanding of chiral symmetry breaking through Random Matrix Theory. 

The purpose of these very elementary lectures was to acquaint students with the new 
developments in our understanding of chiral symmetry breaking in QCD. The amount of time 
was insufficient to give a complete review and this is reflected in these notes which have focused 
on only a few of the many interesting aspects. That chiral Random Matrix Theory can say 
something exact about a complicated quantum field theory such as QCD with light quarks sounds 
like a wild idea, and it was initially met by a lot of skepticism in the lattice QCD community. All 
of this skepticism turned out to be unwarranted. There is now complete understanding of how 
and why chiral Random Matrix Theory represents an exact limit of the light quark partition 
function of QCD in a specific finite-volume regime. A precise mapping can be made between 
observables computed in QCD in that regime, the standard chiral Lagrangian of light-quark 
QCD and chiral Random Matrix Theory. This is a fantastic achievement, a most surprising 
and deep relation between a quantum field theory as complicated as QCD and a sequence of 
universal phenomena that follow from the spontaneous breaking of chiral symmetry alone. One 



example suffices to illustrate the enormous progress that has been achieved. Let us consider the 
eigenvalues of the QCD Dirac operator D: 

Dlpn = X n 1p n (7) 

To make the discussing well-defined, let us consider this eigenvalue problem in a theory with 
a finite ultraviolet cut-off A and a finite infrared cut-off L. As will be discussed further down 
in these lectures, spontaneous breaking of chiral symmetry requires that the smallest non-zero 
eigenvalues A n accumulate towards zero at a rate proportional to 1/L 4 in four dimensions. This is 
thus a simple test: do they or do they not? Before the developments described here, the strongest 
statement that could be made on basic principles was a proof that the A n 's accumulate at least as 
fast as 1/L [TJ. But this is a very weak condition. Even in a free theory, eigenvalues accumulate 
as 1/L. All the theorem says, then, is that Dirac eigenvalues in QCD do not accumulate with 
a rate that is slower than that of free fermions, i.e., as if there were no gauge interactions 
whatsoever. Can one not do better? With the advent of chiral Random Matrix Theory this 
problem has been solved, and it is now known that in QCD with two light flavors (the u and d 
quarks), Dirac operator eigenvalues do accumulate towards the origin at a much faster rate that 
goes exactly as 1/L . Not only that, the precise probability distributions of single individual 
Dirac eigenvalues, Aj, i = 1, 2, 3, . . . ordered according to their distance from the origin can be 
computed exactly. These distributions, as well as all spectral correlation functions of these 
smallest Dirac operator eigenvalues follow universal scaling laws whose exact analytical forms 
are known. 

2. Chiral symmetry in QCD 

If we look at the Lagrangian density of QCD with massless quarks, we are struck by the fact that 
left-handed and right-handed field decouple. If we had exactly Nf massless quarks in QCD, the 
global symmetry group would be U(Nt) x U{Nf). Of this a subgroup is U(l) baryon number, 
which remains unbroken in QCD. The singlet chiral U{1) symmetry is apparent only: it is 
broken by the chiral anomaly, and is therefore not a symmetry of the quantum theory. What 
remain are two independent flavor rotations SU(Nf) for the left and right handed quarks. If 
this symmetry is broken as expected, SU(Nf)L x SU(Nf)r — > SU(Nf) due to the formation 
of a flavor- independent chiral condensate {tpip}, we have Nf — 1 broken generators and hence, 
by Goldstone's Theorem, Nj — 1 massless bosons. This is a profound statement! How this 
slowly came to be realized before quarks had even been introduced is a story in itself. One 
consequence of this is that a low-energy representation of the QCD partition function had been 
inferred long before QCD was considered. In the next subsection we will give a lightning review 
of the moderne viewpoint on this. 

2.1. The Chiral Lagrangian 

Because of Goldstone's Theorem, the low-energy degrees of freedom of massless QCD are those 
of the Nambu-Goldstone bosons. In reality, even the u and d quarks are not exactly massless. 
But on the QCD scale of ~ 1 Gev they look nearly massless; their masses are on the order a few 
MeV. At some point it must be possible to treat the quark masses as small perturbations on top 
of a theory of genuine Goldstone bosons. The chiral Lagrangian does exactly this. Consider first 
a truly massless A^-flavor theory where chiral symmetry is spontaneously broken as discussed 
above. To describe the coset of chiral symmetry breaking it is convenient to introduce an unusual 
non-linear field representation in terms of group elements U(x) = ex.p[i\/2&(x) / F], where F is 
the pion decay constant and &(x) = X a (f> a /V2 represents the collection of Nj — 1 pion fields. A 
field theory based on this should have the following ingredients: (1) kinetic energy terms for the 
pions with canonical normalization, (2) vanishing interactions at zero energy (as follows from 



Goldstone's Theorem) and (3) invariance under the chiral rotations. In the absence of any other 
principle that could fix terms in the Lagrangian the only solution will be to write down all terms 
that comply with the two conditions above. This gives us an endless series of invariants: 

F 2 

C = —Tr[d^d»U] + Li(Tr[d M f7* d^U}) 2 + L 2 Tr[d tM U^d u U]Tr[d v tfd»U] + ... (8) 

with an infinite series of couplings (we display only the first two, L\ and L2 here). Values of all 
these couplings are left totally unspecified. They must depend on details of the dynamics (which 
gauge group is responsible for the symmetry breaking and so on). This is the chiral Lagrangian. 

This cannot be the full story because we also want to treat quark masses, at least in an as yet 
not totally specified perturbation to the above Lagrangian. There is a beautifully simple way 
to do this, which goes under the name of the spurion technique. The chiral Lagrangian ([8]) is 
by construction invariant under global transformations U — > UrUU^. A mass term in the QCD 
Lagrangian, 

tpLmipR + tpRm^ipL (9) 

is of course not invariant under tpi — > Ul^Li^r — > Uri^r. But if the mass (matrix) were a 
field transforming like m — > UlitiU r then the QCD Lagrangian with this mass term would be 
invariant. Then we know what the corresponding terms must look like in the chiral Lagrangian. 
The leading term will be 

£ = T,Tv[mU + m ] U ] ] (10) 

and there will of course be higher order terms involving higher powers of m and derivatives. 
Such a chiral Lagrangian will be invariant under the combined transformation 

m U L mU R , U -»■ UrUU ] l . (11) 

A Lagrangian constructed in this way will thus explicitly break chiral symmetries in precisely 
the same way as the underlying theory of QCD. 

2.2. Different counting schemes 

So far we have defined the chiral Lagrangian by the totally of all terms that transform properly 
under chiral transformations. This looks rather hopeless, since there are infinitely many terms. 
Does such a theory have any predictive power at all? Weinberg [8j showed that it indeed docs. 
First, since all possible terms are already included in the Lagrangian, it is by construction 
renormalizable since quantum corrections cannot produce terms not already present in the 
tree-level Lagrangian. This goes counter to everybody's intuition about renormalizable field 
theories, but it is correct. The chiral Lagrangian is an effective field theory and as such perfectly 
renormalizable. The theory, however, is not well-defined at all energy scales, and this is what 
we should consider next. The trouble is the appearance of higher and higher derivatives in the 
Lagrangian. Normally, this indicates that theory will not be a truly local theory and this may 
hold here as well. Certainly, the more derivatives, the more there will be contributions from the 
high energy scale. But the theory is not even supposed to hold at the high energy scale - it is 
precisely a low-energy effective theory. The way to deal with this was explained by Weinberg: 
consider an ordering of the Lagrangian where one expands in the number of derivatives. If this 
is to make sense, the Lagrangian must necessarily be viewed as a theory with a cut-off. In QCD, 
one would expect the cut-off A to be around 1 GeV; in detail this shows up in the form of the 
combination 4ttF, which indeed is of that order. Dependence on this cut-off can be removed for 
observables below this energy scale. 

The expansion in terms of the number of derivatives can be made precise, see e.g. ref. [9]. 
A counting is introduced which makes it systematic. Since it is in terms of derivatives (or 



momenta), this is referred to as the p-expansion. The pion mass will also enter explicitly in the 
chiral Lagrangian, due to non-vanishing quark masses that explicitly break chiral symmetry in 
QCD. So an ordering with respect to m 2 needs to be introduced as well. The natural expectation 
works: one can systematicall treat as being of order momentum p. Then chiral perturbation 
has only one expansion parameter, p. The scale is given by the expected ultraviolet cut-off 
4irF so that the dimensionless expansion parameter is ~ p/(4irF), This expansion works quite 
well phenomenologically, at least in the sector of the two lightest quarks (and we shall only be 
concerned with this light quark sector here). 

The p-expansion breaks down eventually if one takes the massless limit at finite four- volume 
V. This is most easily seen by considering the pion propagator 

A(p 2 ) = 1 . 1 - . (12) 
y ' V p 2 + ml v ; 

For simplicity, let us take the finite volume V to be a symmetric four-torus of length L. The 
smallest (quantized) momentum is thus of order 1/L. In a massless theory this makes the 
propagator A(p 2 ) vanish for large L like L~ 2 . But for the momentum zero mode the cut-off 
is given entirely given by the mass. Clearly, for m n ^> 1/L the propagator is still protected 
by the mass, but as m n L ~ 1 or much smaller we enter a new regime. From this point and 
onward the usual perturbative expansion cannot make sense in the zero-mode sector. Is this of 
importance? Can we not just ignore the zero mode? The answer is unfortunately no. Although 
it concerns only one mode, this one mode potentially overwhelms all other contributions from the 
perturbative expansion. The resolution of this problem is to treat the zero-mode sector exactly 
in a sense that will be more clear below |10t [TT] . Gasser and Leutwyler developed a modified 
perturbation theory that includes this feature, while retaining all other properties of the usual 
pertubative expansion of chiral Lagrangians. This has become know as an e-expansion (not to 
be confused with the expansion of similar name that expands in dimensionality 4 — d). The 
non-zero momentum modes still retain their usual pertubative expansion since the propagator 
for these modes will go at least as L~ 2 , even in the chiral limit. 

A systematic chiral counting for this so-called e-regime is as follows Let e ~ 1/L. Then 

~ e 2 , p ~ e . (13) 



It is more instructive to think of this in terms of the microscopic degrees of freedom. Let us 
restrict ourselves to two degenerate light quarks of mass m. Due to the Gell-Mann-Oakes- 
Renner relation 

F 2 ml = 2mS (14) 

where £ is the chiral condensate, we see that m ~ e 4 . In the e-regime the quark masses scale as 
inverse powers of the volume. In fact, we can think of this as a regime where the extreme chiral 
limit m — > is taken in a way that correlates with the way the four- volume V is sent to infinity. 
The relevant proportionality factor must have dimension 3, and indeed the right way to think 
of it is that the combination mTiV is kept of order unity. Here the condensate £ provides the 
missing constant of proportionality, which is directly related to the fact that the combination 
rmpip is a renormalization group invariant, in fact it is the explicit chiral symmetry breaking 
term in the QCD Lagrangian. 

There is another way in which we can understand the appearance of the scale mTiV . In 
two-flavor QCD, the chiral condensate (tpip) at finite quark mass m and finite volume V is the 
trace of the quark propagator, 



W> = ^T7d m lnZ 




(15) 

where we have used that each quark contributes to {$ip) with the same amount. Use has also 
been made of the fact that every non-vanishing eigenvalue Xj can be matched by it opposite- 
sign counterpart. This follows directly from the spectrum being chiral, {D,7 5 } = 0. The 
above spectral representation can be used to deduce the so-called Banks-Casher relation for the 
infinite-volume chiral condensate: 

E = lim lim 2m f°° dX , P ^ = = vrp(0) (16) 

m->0V->oo Jo X 2 +m 2 

This is a formal expression that should be understood only in terms of a cut-off theory. The 
chiral condensate {ipip) is ill-defined in the ultraviolet, and indeed one expects the spectral 
density p(X) to behave like ~ A 3 (as in a free theory) for very large momenta, and this would 
lead to a quadratic divergence (ipip) ~ mA 2 . A term like ~ mln(A) is also expected. However, 
the important point here is not the behavior at the ultraviolet end of the spectrum, but near 
the origin. Whether chiral symmetry is spontaneously broken or not is determined by whether 
or not the spectral density p(X) vanishes there. If p(0) is to be non-vanishing, the discrete 
eigenvalues Xj must accumulate there at a rate inversely proportional to the volume [12] : 

AA ~ 1/V near A~0 (17) 

so as to yield a finite density at the origin. The inverse of the proportionality factor is precisely 
the condensate S on account of the Banks-Casher relation. We thus find that XT,V will be 
of order unity for the smallest eigenvalues if we are to generate a chiral condensate. As will 
become clear shortly, there is an e-regime of a partially quenched chiral Lagrangian hidden in 
this statement. 

Before defining the partially quenched theory, let us return to the integral representation of 
(ipip). If we know how this condensate as a function of the mass m, can we not invert the relation 
to get p(A)? The trouble is that p(X) is also a function of quark masses m. If it were not, we 
could indeed simply invert the relation. The trick to achieve this is to do partial quenching: 
one introduces new unphysical quarks that are used to produce an intermediate generating 
function, but which are removed at the end. The spectral density will then be independent of 
the masses of these physical quarks, but if we know the analogue of their condensates (ipip) , we 
can invert the integral representation and compute the spectral density. There are two ways 
of introducing such partially quenched quarks: one is by means of a graded structure [13], the 
other by means of replicas |14^ 115]. Both methods work beautifully in perturbation theory, and 
for the same reason: observables computed in perturbation theory will depend polynomially on 
the number of additional fermionic Nf and bosonic Nf, species. If both fermions and bosons 
are kept, the anticommutation sign from fermions just cancel the corresponding contribution 
from the fermions. So an equal number of degenerate fermions and bosons do not contribute 
to the partition function, a statement that looks entirely trivial if phrased formally in terms of 
functional determinants. In the replica method, one keeps only the fermions (or, equivalently, 



only the bosons) and instead removes them from the partition function by taking the replica 
limit Nf —7- 0. Since the behavior is polynomial in perturbation theory, it is unambiguous how 
to do so in that setting. 

In the e-regime, integrals have to be done exactly, beyond perturbation theory. Curiously, 
then, both the graded method and the replica method run into technical difficulties. In the 
graded method it turns out that the naive extension of the chiral Lagrangian to a graded coset 
fails to yield convergent integrals, a problem that never surfaces in perturbation theory. This 
problem can be overcome by a judicious choice of bosonic integration [IB] : The graded integration 
domain must be a combination of a compact and non-compact regions. In the replica framework 
the difficulty arises from the precise specification of how to analytically continue in the integer 
Nf. In perturbation theory this is trivial (because dependence on Nf is polynomial only), but 
beyond perturbation theory one needs more structure as guidance - knowing a function at all 
integers is of course not sufficient for defining an unambiguous analytic continuation into the 
real line. Fortunately, the required structure is precisely present in the cases of interest [TT] . 
in fact there is a remarkable connection to the theory of integrable systems and the so-called 
Toda lattice hierarchy. It is beyond the scope of these simple lectures to go into that fascinating 
subject, to which we here only refer to the literature [18J. It was based on this structure, and 
its relation to the graded analog, that the partially quenched partition function of the chiral 
Lagrangian was first computed in the e-regime \17\ [T9] . 



3. What is Random Matrix Theory? 

What is a "random matrix"? It is perhaps not obvious, but with the conventional definition we 
must fix its size. Let us first, for simplicity, consider square matrices of size N x N. Elements 
in this matrix can now be chosen at random, according to our chosen distribution. But this 
cannot be the whole story, because some matrices have particular symmetries and particular 
transformation properties under conjugation. Indeed, a sensible definition of Random Matrix 
Theory is based on a division into certain classes of random matrices, a classification that dates 
back to Dyson |20j . The three main classes are labelled by an index /?, which can take values 1,2 
or 4. Here we will focus the /3 = 2 class where one picks the random matrices to be Hermitian, 
and thus having real eigenvalues. 

To produce a Random Matrix Theory one does as in statistical mechanics and sums over an 
ensemble of a given set of random matrices M. This produces a partition function, 

Z = J dM exp [-iVTr(M 2 )] , (18) 

where, for simplicity, we have here restricted ourselves to a Gaussian distribution. We can 
associate with this a probability distribution 

P(M) = -^exp \-NTr(M 2 )} (19) 

and for Hermitian matrices the measure is 

N 

dM = ]]_dM ii ]]_d(RjeM ij )d(JmM ij ) . (20) 

i=l i<j 

We are not interested in the matrices themselves, but rather in their real eigenvalues Aj. It is 
therefore advantageous to go to an eigenvalue representation of the partition function by means 
of diagonalization. This is achieved by a unitary transformation: 



M = U ] DU , D = diag(Ai, . . . , Ajy) 



(21) 



with U*U = 1. 

We note that the Hermitian matrix M is described by N 2 coefficients, while the diagonal 
matrix D has N coefficients and U has N 2 — N. So we can indeed make the change of variables 



M 



D,U 



Since the Haar measure on U(N) is left and right invariant, the asociated Jacobian J of this 
change of variables can only depend on D, i.e., J = J(Xi, ■ ■ ■ , Xn) and 



dM = (j{d\}j J(X 1 ,...,\ N )dU . 



To compute J, it suffices to go near the identity in U(N): 

U = l + ie 

where, for unitarity, e' = e. Then 

M = U^DU = D-i[e,D] 



or, in components, 
In other words, 



Mij = XiSij + ieij{Xi — Xj) [no sum on 



i = j : dMa = dXi 



i < j 



d(ReMij) = (Xi - Aj)d(Re(^y)) 
d(lmMij) = (Xi — Xj)d(Im(iEij)) 



so that 



dM = 



/ N \ 

11 dMa Y[ d(ReMij)d(lmM, 

\i=l / i<j 

\i=i i i<j 



up to the group volume of U(N). Note also that 



iV 



Tr(M 2 ) = Tr(tf DUtfDU) = Ti(D 2 ) = ^ A 



i=i 



(22) 

(23) 

(24) 
(25) 



(26) 



(27) 



(28) 



so that also the probability measure can written entirely in terms of eigenvalues and the partition 
is then 



j \i=i / i<j 



(29) 



There is an easy mnemonic to see that this Random Matrix Theory is in the (3 = 2 class since 
it can be read off from the power "2" in 

i<j 



This particular combination of eigenvalues can be written in the form of a so-called Vandermonde 
determinant raised to the second power: 



n(A,-A,) 2 = A({A 4 }) 2 



(30) 



i<j 



where 



A({AJ) = det(A') j = l,...n, i 



= 0...n- 1 . 



(31) 



The Random Matrix Theory (|29p is of course well defined for all N. When N is finite, the N 
eigenvalues will jump around in a rather random manner. But the presence of the Vandermonde 
determinant (squared) in front of the Gaussian damping factor has an important consequence: 
level repulsion. The probability of two eigenvalues approaching each other goes to zero. For 
finite N, there could be two ways to get around this: 1) all eigenvalues could spread all over 
the real line or 2) eigenvalues could be confined to a finite interval, but level repulsion would 
force their distributions to lock in an essentially equidistant form. This will clearly depend on 
how we rescale eigenvalues. It turns out that the Gaussian damping factor in the probability 
distribution (|29p is so strong that eigenvalues are forced to lie on a finite interval, and it is thus 
the 2nd realization of eigenvalue distributions we are seeing with such measures. The factor of 
N in the exponent is precisely inserted so as to ensure this compact distribution without further 
rescalings of the eigenvalues. 

The distribution of eigenvalues becomes beautifully simple in the limit of very large N: we 
get the famous Wigner semicircle distribution. 

3.1. Chiral Random Matrix Theory 

We are now ready to consider a Random Matrix Theory that is chiral. The identification of this 
class of theories to be the one relevant for low-energy QCD was first made by Verbaarschot in 
a remarkable series of papers about 15 years ago |21l I22j . The idea was to base the Random 
Matrix Theory on transformation properties of the Dirac operator D in QCD. In particular, since 
D anticommutes with 7 s , {D,j§} = 0, a random matrix M is introduced that has a similar 
structure. Moreover, gauge field configurations can be assigned to different classes depending 
on their topology, which is classified according to the invariant 



where F^ u = {l/2)e llua ^F^ 1 ' is the dual QCD field strength tensor. In the QCD partition function 
one sums over all these different classes that have v = 0, ±1, ±2, ±3, . . ., but it can be useful 
to consider also sectors of fixed v. A well known mathematical theorem ensures that in such 
sectors there are exact zero modes of the Dirac operator, all of definite chirality. The number of 
positive minus negative chirality zero modes is precisely the index v. Leutwyler and Smilga [12J 
suggested to study properties of the Dirac operator in such sectors of fixed topology by means 
of the relationship to the chiral Lagrangian. This gives an exciting new opportunity for cross- 
checking how low-energy observables in the full theory (QCD) are in one-to-one correspondence 
with observables in the effective large-distance theory of the chiral Lagrangian. 

Verbaarschot and Shuryak [21] introduced the following chiral Random Matrix Theory: 




(32) 




(33) 



where V(M 2 ) will be discussed below, and where 




is a {2N + u) x (2N + v) block Hermitian matrix composed of the complex N x (n+v) rectangular 
matrix W. The integral is over the elements of this matrix. We note that both measure and 
integrand is then invariant under 



W -> V^WU , U 6 U(N) , V e ?7(iV + v) (34) 

and for this reason this Random Matrix Theory is called chiral unitary. But why is it 'chiral'? 
Consider the (2N + v) x (2N + v) matrix 



7 5 
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This matrix clearly anticommutes with the matrix M, 

{M, 7 5 } = , (35) 

which implies that all non-zero eigenvalues of M come paired ±Aj. So the matrix M shares 
the U(l) chirality property of the full Dirac opertator of QCD. Moreover, due to its rectangular 
block decomposition, there are v zero modes. The spectrum of M therefore has two parts: there 
are 2N non-zero eigenvalues that are chirally paired and v zero modes. 

Actually, the chiral Random Matrix Theory (|33p was originally introduced with just a 
Gaussian measure, but it was believed from the beginning that the so-called 'microscopic' results 
that could be derived from it would be universal, i.e. to a very large extent independent of the 
chosen potential. Indeed, it was proven in ref. |23j that all results to be discussed below hold 
for an arbitrary potential 

V{M 2 ) = J2^rM 2k (36) 
k>i k 

with, essentially, only one single constraint: the spectrum of eigenvalues of M must have support 
at the origin: p(0) 7^ 0. This is in beautiful accord with the intuition, to be be made more precise 
below, that a non-vanishing chRMT spectral density at the origin is required for this theory to 
describe aspects of chiral symmetry breaking - in heuristic analogy with the Banks-Casher 
relation in QCD. The universality and robustness of this result is crucial for the understanding 
of why chiral Random Matrix Theory can describe chiral symmetry breaking in a full-fledged 
quantum field theory such as QCD. 

Interestingly, one can also consider a chiral Random Matrix Theory where the spectral density 
at the origin p(0) is carefully tuned so as to produce a zero: p(0) = [24j . One might hope that 
this could represent chiral symmetry restoration, one way or another. This, however, has never 
been established. There are still universality classes in this case, one for each order of zero of 
the spectral density. Each class can be reached, successively, by tuning more and more of the 
parameters in (|36p . This is an example of universality in what can be dubbed multicriticality 
(while the more general condition p(0) 7^ corresponds to the non-critical case). A pity that 
there seems to be no known instances of quantum field theory where these more sophisticated 
universality classes can be of relevance! In any case, the mere existence of these extended classes 
shows that the universality domain of the non-critical ensemble is restricted and by no means 
trivially implied for any chiral Random Matrix Theory potential V(M 2 ). 

Very similar to the way we could introduce an eigenvalue representation (|29[) for the non-chiral 
(Gaussian) ensemble, one can also find an eigenvalue representation 

N , 
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where the zf's are the real (and positive) non- vanishing eigenvalues of W^W, and where in the 
last line we have just transformed to more convenient variables Aj. We have also used the fact 
that diagonalization of W^W automatically gurantees factorization of the potential part in the 
manner shown. It is interesting to note how the dependence on Nf and v has merged into the 
combination Nf + u only. This is a simplifying property that holds only in the massless case. The 
case of massive quarks will be discussed below. The square of the Vandermonde determinant 
has appeared just as in the non-chiral case. Also here it will result in an eigenvalue repulsion. 

There are standard methods for computing the spectral correlation functions of a Random 
Matrix Theory with an eigenvalue representation such as (|37p . A convenient method is based 
on orthogonal polynomials, using a beautiful formalism that dates back to Dyson. Limitations 
of space does not permit a detailed discussion of this, and here we will restrict ourselves to 
explaining the main results that follow from that kind of analysis. 

Let us from now on work entirely in rescaled variables £ = 2ttN p{Q)\. For fixed Q this implies 
smaller and smaller eigenvalues - those near the origin. A central object is the kernel i^(Ci 5 (2)- 
If we know this kernel we know all spectral correlation functions 

p(( 1 ,...,(k) = detK{( i ,Q (38) 
and in particular the spectral density itself, 

p(() = K(C,Q • (39) 
The kernel appropriate for QCD is the Bessel kernel derived by Verbaarschot and Zahed [22J, 

-Ci-WCiW Cg) - C2 J a (Ci)^+i(C2) 

d 2 - C2 2 



*(Ci, Gi) = y/C& * lua + 1KSLJUa ™ ^^±^1 f (40) 



where a = Nf + v (we consider only positive v here). Despite appearance, this is regular at 
coincident points where it gives the microscopic spectral density 

P(0 = ~ (Ja(() 2 ~ J a+ l(C)Ja-l(C)) ■ (41) 

This is the most simple and remarkable result of the chiral Random Matrix Theory approach. 
The claim is that this gives the microscopic spectral density of QCD with Nf massless quarks 
in the finite-volume regime we consider. Below we shall see why this is so, and how to include 
massive quarks as well. Amazingly, if we boldly set Nf = this should give the spectral 
density of the Dirac operator in pure Yang-Mills theory! When we see the connection with the 
partially quenched chiral Lagrangian this makes it clear why this should be possible (with some 
caveats). Another immediate check on the above spectral density comes from the fact that it 
reproduces exact spectral sum rules derived by Leutwyler and Smilga [12] on the basis of the 
chiral Lagrangian alone. 



3.2. Spectral properties of the Dirac operator 

The pioneering paper that showed how these detailed predictions for the spectral density of 
the Dirac operator can be tested in lattice gauge theory simulation was ref. [25]. Quenched 
staggered fermions with gauge group SU{2) was used there. Because the chiral symmetries of 
staggered fermions away from the continuum limit are different from those of continuum fermions 



when the fermions carry real or pseudo-real representations of the gauge group (pesudo-real in 
the case of fundamental fermions and gauge group SU(2)), the detailed predictions were not 
those of eq. (|4ip . But also for this case the analytical prediction had been worked out [26] and 
although it differs in detail, it shares many of the properties with that of ()4ip . The analytical 
expressions for continuum fermions based on this gauge group and fermions in the fundamental 
representation had been worked out as well |27| . We shall not go into great detail with these 
different predictions, but it is an amazing fact, first recognized by Verbaarschot [28J, that there is 
a one-to-one correspondence between the Dyson classification of (chiral) Random Matrix Theory 
and the three chiral symmetry breaking patterns in vector-like gauge theories: 

• The fermion representation r is pseudo-real: Chiral symmetries are enhanced from 
SU(Nf) x SU(N f ) to SU(2N f ), and the expected symmetry breaking pattern is SU(2N f ) -> 
Sp(2N f ). 

• The fermion representation r is complex: The expected symmetry breaking pattern is 
SU(N f ) x SU(Nf) -> SU(N f ). 

• The fermion representation r is real: Chiral symmetries are again enhanced to SU(2Nf), 
and the expected chiral symmetry breaking pattern is SU(2Nf) —> SO(2Nf). 

That there is indeed a one-to-one matching with the three Dyson classes, and that the chiral 
symmetry breaking patterns really do follow the above classification, has been tested on the 
lattice in great detail for a variety of different fermion representations and it works in all cases 
[29]. 

There is one further general peculiarity of simulations with staggered fermions on the lattice 
which in the beginning caused some confusion, but which really is just one further confirmation 
of the above classification scheme. This arises from the fact that staggered fermions away from 
the continuum limit have an additional U(l) factor in the coset of chiral symmetry breaking 
(this holds for all three classes). As we will show below, this means that all results for staggered 
fermions away from the continuum are equivalent to those we would get if we had projected on 
the v = sector. Any attempt at a further projection on sectors with v ^ leads to the same 
partition function. So for staggered fermions away from the continuum it will look as if only 
the sector of v = appears. This has been confirmed by a high-statistics study where gauge 
field configurations were even selected according the their (approximative) topological indices: 
all spectral data of the staggered Dirac operator look as if they have been obtained from the 
v = sector [30J. This should of course change as the continuum limit is reached. By using 
highly improved actions as well as improved Dirac operators that smear both configurations 
and observables it has been checked that predictions of continuum fermions do match those of 
staggered fermions eventually [31] . 

While most initial studies of the spectral properties of the Dirac operator were restricted 
to gauge group SU(2), results were soon available for the QCD gauge group of SU(3) [32] . 
nicely confirming the prediction (|4ip . A far more exciting possibility was to check all of these 
predictions with 'chirally good' fermions based on the overlap operator. This was first done in 
ref. [33]. Here the effect of topology can be clearly seen, and this works well even for approximate 
realizations of the overlap operator [M] . A high-statistics study with several eigenvalues included 
has been made in ref. [36], showing clear agreement with the predictions. It is also interesting 
to see how, with an infinite number of colors, chiral symmetry can break even at finite four- 
volume. Results still fall right on the chiral Random Matrix Theory predictions, now as a result 
of scaling in the number of colors N c [37] . Here one sees clearly how color and volume acquire 
complementary roles at large- N c , the corresponding scaling variable being \J}N C V. 

So far our discussion has been restricted to the massless limit, which is easily taken in lattice 
gauge theory if one restrict oneself to the quenched theory. An interesting question is how to 
include quark masses in all of these preditions. If we consider the chiral Random Matrix Theory 



()33|) . it is clear that masses m must be rescaled at the same rate as eigenvalues Aj. Let us 
re-emphasize that m in (|33p is simply a dimensionless parameter in the matrix integral. It is 
not a quark mass. However, just as the universal predictions for rescaled eigenvalues Np(0)X 
correspond to volume-rescaled eigenvalues VT,X on the QCD side, so Np(0)m will correspond 
to the scaling variable VEm in QCD. The precise predictions were worked out in ref. [39J. The 
clue to the solution was there found in the theory of orthogonal polynomials, which is known to 
solve the spectral problem of eigenvalue representations like those of (|33|) . Intuitively, it comes 
about as follows: if we know a set of orthogonal polynomials corresponding to (|33p without 
the determinant in the measure, we can construct the set of otrhogonal polynomials with one 
determinant in the measure. Basically, one expands the new set of orthogonal polynomials in 
the basis of the old set. The coefficient in this expansion can be fixed by what is known as 
Christoffeks Theorem. This solves the Nf = 1 case in terms of the quenched case. To illustrate, 
one finds the spectral density in sector of topological [39] charge v = 0, 

p(C.m) = I {MC? + MC?) - < *(0fl Wff'O + W WO] 

2 v ' (C + /U 2 )W) 

In this way of writing it, we nicely illustrate an essential feature: 

}mip((, t i)=p(0 (43) 

in other words: when we take the (rescaled) quark mass to infinity we recover the quenched 
spectral density. This is as we would intuitively expect, since the quark should then decouple 
and the eigenvalues should distribute themselves according to the quenched distributions. 

All of this was generalized in ref. [39] to any number of flavors, and the universality proof 
was there extended to this generalized case as well. Moreover, the explicit expressions for the 
massive kernels were provided there. From these, all spectral correlation functions with massive 
quarks follow. However, the expressions can become unwieldy, and there is actually a much 
more compact way of understanding these results. It turns out that, quite surprisingly, one can 
expression all microscopic spectral correlators compactly in terms of finite volume QCD effective 
partition functions. The crucial ingredient is a powerful identity derived by Zinn-Justin [3D] in 
the case of the ordinary unitary ensemble, and which is readily generalized to the chiral case. 
Let us consider the v = case in detail JUJ. The kernel can then be expressed as 
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The last integral is over (N—l) eigenvalues only. However, in the large- N limit we shall consider 
below, this distinction can be ignored. Thus, in the large- N limit we have 



K N (z,z';mi,...,m Nf ) = e~% ^ 2 )+V(z")) y ^pj^ /^ 2 + m j)(z' 2 + mj) 
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Z v (mi, . . . ,m Nf ) 



where the matrix model partition function in the numerator is evaluated for a theory 
corresponding to (Nf + 2) fermions, of which two have imaginary mass. By means of the usual 



factorization property, all higher n-point spectral correlation functions are then also explicitly 
expressed in terms of the two matrix model partition functions Z^ 1 ^ and zl Nf+2 \ The spectral 
density corresponds to the two additional (imaginary) masses being equal: 



p( N f\z;mi,...,m Nf ) = lim K N (z, z; mi, . . . , m Nf ) • (46) 

J iV— KX) 1 

All of this is exact, and at the level of unrescaled variables. We can now take the microscopic 
scaling limit in which ( = zN2irp(0) and pi = miN2irp(0) are kept fixed as TV— >oo. In this limit 
the pre-factor exp[— (N/2)(V(z ) + V{z' 2 ))] becomes replaced by unity. Identifying £ = 2-7r/)(0), 
this is the limit in which we can compare with the finite-volume partition function of QCD. 
Proceeding in this way, we get, with C being an overall constant, 

(N f +2) ( ... 

P ^\Cwi,... t m t ) = c\c\l[(e + p}) Zv ^•••>/^XX) (47) 

/ Z v (p 1 , . . .,p Nf ) 



and all n-point correlation functions are given by one magical Master Formula [41] 

p ( N f>») (Ci, - - . , Cr»; Mil - - - j mjv/) = det K(( a , C&; ■ ■ • , i*N f ) • (48) 

a,b 

This generalizes to any v [42j and in fact leads to some surprising identities among the partition 
functions involved, identities that express the connection to an associated integrable hierarchy. 
The advantage is that the finite- volume partition functions are known in exact analytical forms 
|43| . Results have also been generalized to the other chiral symmetry breaking classes [441 143] . 

One immediate check on these massive spectral densities is that they reproduce exact massive 
spectral sum rules [46] just as the massless spectral densities reproduce the massless spectral 
sum rules of Leutwyler and Smilga. 

Amazingly, the same procedure generalizes to the computation of individual eigenvalue 
distributions. Such a notion makes clear sense in the chiral Random Matrix Theory framework 
since there is a natural starting point from where to start counting eigenvalues (the origin). 
After the lattice gauge theory community had accepted that chiral Random Matrix Theory 
gives exact results for the spectral density of the Dirac operator, and all higher-point spectral 
correlation functions, there remained a curious belief that even the very notion of ordered Dirac 
operator eigenvalues and their precise distributions could not be given field theoretic meaning. 
This is of course not so: ordered Dirac operator eigenvalues have precisely as much meaning in 
the field theory framework as the spectral density and all spectral correlation functions. This 
will become clear below. But as first step towards this, we should describe how individual 
eigenvalue distributions of the chiral Random Matrix Theory can be given explicitly in terms 
of finite volume partition functions for first the smallest Dirac operator eigenvalue [47] and 
then the general distribution for the kth Dirac operator eigenvalue [48J, as counted from the 
origin. Let us denote the joint probability distribution of the smallest k ordered eigenvalues by 
w fc(Ci) • • • j Ofc-li Cfc)> where £1 < . . . < Ck- By manipulations very similar to the ones used above, 
one can show that it can be written in terms of finite volume partition functions 
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where the "masses" y Ck ~ Cf are doubly degenerate (and 1 < i < k — 1). Similarly, the set of 
"masses" {Cfc} is v-iold degenerate. 

A very high-accurate Monte Carlo study of these individual eigenvalue distributions has 
been done in ref. [49] . for both the quenched case and for theories with dynamical finite-mass 
fermions. The agreement is quite spectacular, with parameter-free fits that lie right on the 
analytical predictions. Of course, if one for some reason would fear these predictions, one could 
also just look at the chiral condensate which follows a very precise scaling law based on these 
analytical predictions for the chiral condensate [50|, l5"T] . 

One important question is how individual eigenvalue distributions, ordered in the way we 
have discussed above, can be understood and derived from the quantum field theory point of 
view. How do we even define such distributions in quantum field theory? It turns out that 
there is a simple answer |53j . For simplicity, let us here restrict ourselves to the case of just the 
smallest Dirac operator eigenvalue, the distributions of the larger eigenvalues follow analogously 
|53j . Here one can show that the so-called gap probability has a convergent expansion of the 
form 

Eo(s) = 1 - fdX 1 p(Xi) + \ f dAidA 2 />2(Ai, A 2 ) + • • • (50) 
Jo * Jo 

involving higher and higher spectral correlation functions. All of these spectral correlation 
functions can be computed in quantum field theory by means of suitable sources. The 
distribution of the smallest eigenvaluie now follows by differentiation: 

d 

pi(s) = - g- s E (s) • (51) 

For the kth. eigenvalue this follows by a simple iterative scheme that involves increasingly higher 
order correlation function as the first term in the expansion [53]. But the convergence is fast in 
the appropriate region of A, and this scheme is in fact a quite practical way of computing these 
distributions if one is only interested in high numerical accuracy. 

Another issue concerns gauge field topology. It is clear from the chiral Random Matrix 
Theory construction that topology is intimately tied together with whole construction. While a 
summation over topology is straightforward in QCD, there seems to be no way to peform such a 
summation in the context of chiral Random Matrix Theory. Every sector corresponds to matrices 
of different sizes. However, all physical observables can be decomposed into their counterparts 
in distinct topological sectors, and a summation over topology is therefore indirectly possible 
through [52J 

(O) = \Y. Z ^°)v (52) 

V 

where Z is the full partition function and subscript v indicates when the average is taken in a 
sector of fixed topological charge. Note how averages are weighted by the partition function in 
fixed sectors. Using this relationship and the Master Formula derived above, it is possible to 
perform the sum analytically and derive an analogous formula for the spectral density of the full 
theory |52j . It is also particularly illuminating to see analytically how the triviality of results 
in the Nt = 1 theory are recovered from the highly non-trivial predictions in sectors of fixed 
topological charge. 

3.3. Matching chiral Random Matrix Theory with the chiral Lagrangian 

The Master Formula (|48p gives a strong indication that all these chiral Random Matrix Theory 
results can be derived entirely from chiral Lagrangian framework. This is indeed the case |54|. 1 16| . 
The framework is partially quenched chiral perturbation theory, which we have already discussed 
above. Through it, one can derive a one-to-one correspondence between the spectrum derived 



from the chiral Lagrangian and the spectrum derived from chiral Random matrix Theory. The 
final steps showing that this holds for all spectral correlation functions were taken by Akemann 
and Basile in ref. |55j . 

What has been achieved is thus a precise and detailed understanding of how chiral Random 
Matrix Theory can yield exact results. It provides an exact representation of the leading-order 
chiral Lagrangian in the extreme finite-volume e-regime. The proof goes through the partially 
quenched theory which here serves two purposes at one time: (1) it provides predictions for 
genuine partially quenched lattice gauge theories in this finite-volume regime, and (2) it serves 
as a generating function of spectral correlation functions. In chiral Random Matrix Theory one 
does often not make use of a partially quenched (graded) version of the theory because other 
simpler technqiues (such as the one based on orthogonal polynomials) are available. It is precisely 
by going through the graded version that the full equivalence has finally been established for all 
correlation functions. 

An interesting question is whether one can extract information about the other low-energy 
constants F as well as all higher-order coupling constants of the chiral Lagrangian by means of 
spectral properties of the Dirac operator. In principle, this should be possible since the influence 
of these couplings on the distributions of the Dirac operator eigenvalues can computed. By a 
careful comparison to numerical lattice data, one should be able to extract these other low- 
energy parameters of QCD by considering only a few of the lowest Dirac operator eigenvalue 
- a remarkable thought. However, in practice it is not that easy. The simplest example is the 
leading correction to S, which we already discussed above. In principle, by considering a variety 
of different volumes V, and perhaps different geometries as well, one can extract the pion decay 
constant F from this. 

Other chiral Random Matrix Theory approaches, aimed specifically at determining F, have 
also been invented. One consists in introducing a new source in the chiral Lagrangian that 
couples directly to F to leading order in the chiral counting, while retaining anti-Hermiticity 
of the Dirac operator in that context. Although it is to be considered only as a source, it 
actually has a physical interpretation as imaginary isospin potential. With an external isospin 
potential the QCD measure is free of the sign problem that plagues QCD with baryon chemical 
potential. Making the potential purely imaginary ensures that the Dirac operator eigenvalues 
with this type of chemical potential remain anti-Hermitian. This approach was first developed 
in the chiral Lagrangian framework |56| . but inspired by a construction due to Osborn |57| who 
introduced a two-matrix chiral theory do describe baryon chemical potential, a chiral two-matrix 
theory was set up to describe imaginary isospin chemical potential [58j . Indeed, results from 
this theory, which we will describe next, agree precisely with those obtained from the chiral 
Lagrangian [SB]. Because all the powerful machinery of the more usual chiral Random Matrix 
Theory is preserved, this formulation allows for an explicit calculation of all spectral correlation 
function. 

The first step consists in considering the eigenvalues of two different Dirac operators in QCD, 
Drf™ = MA) + iwo}^ = iA< n) V?° 

d 2 # ee [D(A)+^ 270 h4 n) = i4 n) v£°- (53) 

Imaginary isospin potential corresponds to fi ee pL\ = —fj^. Such a source couples directly to F 
is the chiral Lagrangian for which the leading-order terms in a suitably defined e-regime read: 

Z (*t) = f dU{detU y e ^Tr[U,B][UKB]+^VTr(M^+MUi) (M) 
Ju(N f ) 

where the matrix B =diag({/xi}, {^2}) contains the two chemical potentials introduced above, 
and A4 = diag(mi, . . . , uin,) is the mass matrix. Chemical potential jjL\ has been assigned to 
N\ of the quarks and chemical potential /J>2 to the remaining Nf — N\ = N2 quarks. 



A two-matrix chiral theory that corresponds to this is 
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The matrices and ^ are complex rectangular matrices of size N x (N + v) in the same manner 
as in the usual chiral ensemble. 

Without going into the technical details, an eigenvalue representation can also be found for 
this theory. Let us define 

ci = (l + /^)/<5 2 , c 2 = (l + ^)/5 2 
d = (l + /^ 2 )/<5 2 , 1-r = d 2 /(cic 2 ) 

J = fa- fj,i , (56) 

Then, 
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The two sets of eigenvalues and j/j are obtained after diagonalizing 

$i = $ + 

$ 2 = $ + /i 2 ^, (58) 

and due to this redefinition the original matrices now become coupled in the exponent. This looks 
horribly complicated, but it actually has some recognizable structure. First of all, the parameter 
5 is what measures the strength of the external source, the imaginary chemical potential. Next, 
by a stroke of luck, a formalism of bi-orthogonal polymials had recently been invented for 
precisely this kind of problem by Eynard and Mehta [59J. Analytical miracles occur, and one 
ends with very simple analytical expressions for all spectral correlation functions. Particularly 
interesting spectral observables are those that vanish if 5 = 0. This leads to a strong signal with 
which the pion decay constant can be measured (since 5 couples to F in the chiral Lagrangian). 
All results derived from the chiral Lagrangian in |56j are reproduced in this Random Matrix 
Theory approach. Moreover, all spectral correlation functions can now be expressed in closed 
analytical forms [58]. Also individual eigenvalue distributions can be computed analytically in 
this two-matrix theory [60]. Detailed comparisons with lattice gauge theory data, for a variety 
of different volumes and including the finite- volume corrections |61] . have been made recently 



3-4- Beyond chiral Random Matrix Theory 

Once one has realized the connection between the leading-order chiral Lagrangian in the e- 
regime and chiral Random Matrix Theory, it becomes obvious how to combine the two. The 
most difficult part of the e-expansion is precisely the non-perturbative contribution from the 
momentum zero modes of the pesudo-Goldstone bosons. We have seen how a partially quenched 



chiral Lagrangian had to be understood in this context just in order to derive properties of the 
spectral density. This is extremely useful for lattice gauge simulations since there one is often 
interested in doing partial quenching as an approximation to a real simulation. For example, 
one might wish to scan more parameter values with the same lattice configurations, or one has 
difficulty simulating with light enough fermions. Here partial quenching can provide additional 
information. 

The first issue concerns pion-loop contributions to the leading-order chiral Lagrangian. This 
was investigated in the original paper by Gasser and Leutwyler [11 j. Considering the contribution 
to the chiral condensate, they found that the Lagrangian remains form invariant except for a 
rescaling of the S-parameter, 
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for a four-dimensional torus of volume V. The coefficient is one of an infinite sequence of shape 
coefficients that appear in the e-expansion [63] . This rescaled low-energy constant £ e / t appears 
naturally in the expansion of any observable in the e-regime. We see that the correction term 
disappears for large volume as 1/y/V, which is of order e 2 , as expected. The expansion thus 
seems to work well. In the fully quenched case this is not the case. There a new scale, set by the 
ry'-mass in physical terms, enters and there is no longer a simple one-parameter expansion. As 
usual, when things can go wrong, they will go wrong. In this case one ends up with a new term 
that grows logarithmically with volume - a "quenched finite volume logarithm" [M]. Although 
such a correction term can be kept under control if one only compares two different volumes 
(at least to that order), it is just one more indication of the fundamental difficulty of quenched 
theories. In a physical theory with dynamical quarks such a logarithm still shows up at higher 
orders, but always suppressed by additional inverse powers of L, and therefore harmless. 

Hansen [65] was the first to explore the e-regime beyond leading order for mesonic correlation 
functions in chiral perturbation theory. This was before Leutwyler and Smilga [12] had shown 
the usefulness of working in sectors of fixed topological charge u, and a re-analysis of Hansen's 
results in such sectors has been done [BB]. The generalization of Hansen's results to the quenched 
case was also done there. What is interesting about mesonic correlation functions in the e-regime 
is their polynomial dependence on the euclidean length scale. In the p-regime correlators have 
the conventional exponential fall-off, but in the e-regime, where one is essentially "inside the 
pion", the correlation function is almost flat. One can understand this very simply. Consider 
the zero-momentum projection of the standard massive pion propagator, 
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where T is the (arbitrarily chosen) euclidean time extent. Apart from the overall factor of 1/m 2 , 
we can Taylor expand this in powers of m 2 to find 
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The finite-order polynomials of this expansion are precisely what in the e-regime replace the 
ordinarily exponentially decaying correlation functions. The pole-term at m n = is the 
contribution from the momentum zero-mode, which indeed blows up in the simple expansion 
above (but which is treated exactly in the e-expansion). This is the way in which the e- 
regime expansion of correlation functions merge smoothly with the standard p-regime correlation 
functions. Such an approach has also been developed further, see, for instance, ref. [67j . 



An interesting question is whether one can combine the e-regime with the p-regime, either 
partially quenched or not. Similarly, one can consider the question of so-called 'mixed actions' 
in this context (partial quenching where the light quark limit is reached by using configurations 
generated with different types of fermions, for example). Ideally, one would like to be able to do 
only one calculation which would be valid in a domain that stretches between the two different 
counting regimes. This is an active area of research at the moment [68} 1691 170] and results are 
currently being compared to large-scale simulations based on overlap fermions |71j . 



4. Conclusions 

There is a fascinating interplay between the chiral Lagrangian and chiral Random Matrix Theory. 
Here we have reviewed some of the exact relations that have been established. The fact that 
chiral Random Matrix Theory provides the leading contribution to the e-expansion of chiral 
perturbation theory has turned out to give an analytically very powerful new view on low-energy 
QCD. 

In lattice gauge theory, where most of the results discussed have very direct applicability, one 
is dealing with a discretized version of the continuum field theory. This induces lattice artifacts 
that are not physical. In a first approximation such discretization errors can be neglected, but 
in a more careful analysis they must be taken into account. In fact, measuring the quantitative 
effects of such lattice artifacts can be done in a systematic way following what is known as the 
Symanzik program of effective lattice field theory. There is then a modified chiral Lagrangian 
that includes discretization effects, but phrased in the continuum language. One can then also 
ask for the effect of discretization errors on the Dirac operator spectrum [72]. Remarkably, chiral 
Random Matrix Theory can be extended to include such effects as well [751 EH- It is particular 
interesting to see how the special case of Nf = 1, where there is no spontaneous breaking of 
chiral symmetries, looks in this context. The full partition function is, to leading order, a trivial 
exponential of the free energy since there is a mass gap even in the chiral limit. But it contains 
now two terms: one proportional to the quark mass m, and one proportional to the leading 
lattice-spacing artifact a 2 . These two terms compete trivially in the full theory for which the 
leading term in the effective partition function simply reads 

Z = exp [mSF - 2W 8 y a 2 j (62) 

where Wg is a new low-energy constant of Wilson fermions in this curious Nf = 1 theory (the 
constant of 2 in front is a convention inherited from the general-iVy theory where one has a 
genuine chiral Lagrangian). If one splits the partition function into sectors of fixed index u, 
there is a non-trivial interplay between the two terms above due to the fixing of the index |75j . 
It is quite amazing how much carries over from the case of spontaneous chiral symmetry breaking 
once the index v has been fixed. For the more physical case of a larger number of dynamical 
quarks also other observables, such as space-time correlation functions, can be computed for 
lattice gauge theory observables in the e-regime with explicit a 2 corections. A variety of different 
scaling regimes can be considered. Such a program is presently being carried out for the case of 
Wilson fermions [76] . This is one of the clearly identifiable directions in which one can envision 
future research within this field. 
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